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OFF-NADIR ANTENNA BIAS CORRECTION USING 
AMAZON RAIN FOREST a° DATA 

I.J. Blrrei’ , E.M. Bracalente , G.J, Dome®, J. Sweet , 

G. Berthold* and R. K. Moore* 

ABST R ACT 

The radar response from the Amazon rain forest has been studied to 
determine the suitability of this region for use as a standard target 
to calibrate a scatterometer like that proposed for the National Oceanic 
Satellite System (NOSS). Backscatterlng observations made by the 
SEASAT-1 Scatterometer System (SASS) show the Amazon rain forest to be 
a homogeneous, azlmuthally-i sotroplc, radar target which Is Insensitive 
to polarization. The variation with angle of Incidence may be adequately 
modeled as “ a6 + b with typical values for the Incidence-angle 

coefficient from 0.07 " 0.15 dB/deg. A small diurnal effect occurs, with 
measurements at sunrise being 0.5 dB - 1 dB higher than the rest of the 
day. 

Max 1 mum- 1 1 kel 1 hood estimation algorithms presented hero permit deter 
mlnatlon of relative bias and true pointing angle for each beam. Specifl 
Implementation of these algorithms for the proposed NOSS scatterometer 
system Is also discussed. 


University of Kansas Center for Research, Inc., Lawrence, Kansas 
4 * 

NASA Langley Research Center, Hampton, Virginia 
“Bell Laboratories, Holmdel, New Jersey 
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1.0 INTRODUCTION 

The National Oceanic Satellite System Scatterometer (NOSS SCATT) was 
designed to determine wind vectors on the ocean's surface by measuring 
radar backscatter from the surface. The SCATT was to be a fan-beam scatter- 
ometer utilizing six dual-polarized (V and H) antennas. 

To achieve maximum wind-vector accuracy careful cross-calibration 
between anten.^as must be performed. This report discusses the development 
of a calibration technique utilizing measurements of the scattering coeffi- 
cient 0 ° from the Amazon rain forest. A brief discussion is presented of 
the method used to discriminate rain-forest signals from others. Haps (used 
to screen Brazil data) generated for this purpose are on a 0.25“ x 0.25“ 
grid for locating rivers and a 0.5“ x 0.5“ grid for Identifying vegetation 
types . 

The suitability of the Amazon rain forest for use as a calibration 
target is shown in Section ^.0. From July to October 1978 numerous back- 
scatter measurements were made by the SEASAT-1 Scatterometer System (SASS). 
These measurements show the rain forest to be a homogeneous, azimuthal ly- 
isotropic radar target which is insensitive to polarization. Variation of 
0° with respect to incidence angle may be adequately modeled with a straight 
line fit of the form: = a0 + b , where 0 is the angle of Incidence 

(deg.). Measured values for the Incidence angle coefficient, a, range from 
0.07 “ 0.15 dB/deg depending on beam and time of day. A small diurnal effect 
was observed with measurements made at sunrise ( 050 O-O 630 ) being consistently 
0.5 " 1.0 dB higher than the rest of the day. 

To cross-cal ibrate the beams on a multibeam scatterometer, two para- 
meters, a relative bias and the true antenna pointing angle, must be deter- 
mined for each beam. Development of two maximum-likelihood-estimation 
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algorlthms to determine these parameters using Amazon rain forest data 
is discussed in Section 5.0. The first algorithm is used to determine 
both relative bias and true pointing angle. A second and considerably 
simpler version is then developed to monitor long-term changes in rela- 
tive bias due to drift In the transmitter power. Specific processing re- 
quired to Implement these algorithms for the SCAlT Is also discussed in 
Section 5.0. 

2.0 BACKGROUND 

During the SEASAT program, the need for In-f light cross-calibration for 
antennas from a multibeam scatterometer became apparent. The originally 
proposed procedure using aircraft underflights proved Inadequate, Based on 
the consistency of o° data over the rain forest obtained during SKYLAB, 

R. K. Moore of The University of Kansas suggested using the Amazon rain for- 
est data as a means of calibration. 

The accompanying map (Figure 1) shows the vast extent of the Amazon 

2 

rain forest, it includes an area of more than 3,000,000 km . Because of 
Its equatorial location, seasonal effects are minimized and vegetation lu 
always present. The region is very flat and uniformly forested except for 
breaks in the canopy over the Amazon River and its principal tributaries. 

The predominant climate is hot and humid, with annual precipitation as high 
as 3000 mm. A slight decrease In rainfall is observed In the large central 
part of the rain forest during the months of October to December. A similar 
decrease takes place in the northern part during the months of June to 
September. In the northwest a more significant decrease in rainfall occurs 
during May and June with the lowest rainfall In July and August. 
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FIGURE 1: Extent of Brazilian Rain Forest 
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During the SKYLAB SL-2 mission vertlwl-^polorlzutlon measurements were 
made over the Amazon rain forest at Incidence angles of 32“ to 36 *, Figure 
2 shows a sample time history taken during one of the SKYLAB passes over the 
eastern part of the forest, The scattering coefficient was stable with a 
maximum deviation of aj^out 0,2 dB. For tho combined set of SKYLAB measure- 
ments the mean value of 0 ° was reported as - 5,9 dB with a standard devia- 
tion of I dB [1,2], This observed consistency prompted the further Investi- 
gation using the SEASAT scatterome; -r. 

3.0 DEVELOPMENT OF HAPS USED TO CLASSIFY BRAZIL TARGETS 

3 . 1 Land-Water Map 

For the early analysis of SEASAT observations of Brazil, discrimination 
was achieved by usino vegetation maps developed during the SKYLAB program [2], 
The 1“ X 1“ grid a jle was thought to be adequate. However, as noted 
In Section 4,2, the presence In the SASS footprint of the Amazon River or 
one of Its principal 'tributaries markedly lowered the radar cross-section. 
Since the SCATT was planned to operate with an even smaller footprint, a 
finer map grid was needed to screen out data contaminated by the presence of 
rivers. 

For NOSS a digitized map was created with a 0.25“ x 0.25“ grid. The 
source map was produced by the Department of Cartography of the Brazilian 
Institute of Geography in 1971. This map of the entire Amazon basin Is 
scaled at approximately 3 inches per degree. The digitized map was created 
by overlaying grids to subdivide the source map Into 0.25“ boxes. Codes 
were then assigned to each box using the following conventions; 




"7- 


0 - flatland only (no significant rivers) 

1 - some small rivers but mostly land 

2 - large rivers through box 

3 “ rough terrain 

The results of this procedure were then encoded Into a computer algorithm. 

For a given latitude and longitude the algorithm will fetch the code of the 

box containing the location. This map was tested during the extensive analy- 
sts of SEASAT data discussed In Section If.O. 

For NOSS It was envisioned to use the computed location of the center 
of the Doppler cell for classification. A more careful screening could be 
performed by using the latitudes and longitudes of the corner points for 
a particular Doppler cell. Testing to be sure that all four corners showed 
a "land-only" condition would be a check against overlapping part of a 
rlver-f 1 lied box. 

3.2 Vegetation Map 

To further classify Brazilian data a vegetation map for the Amazon basin 
was produced with a 0.5“ x 0.5“ grid. Codes were assigned to each box accord- 
ing to the conventions listed In Table 1. Six distinct types of vegetation 
are found In the Amazon basin. Types 1 and 2 are the forests which are suit- 
able for use as a standard target. Combinations of types were only designated 
If the dominant vegetation was either type 1 or type 2. Otherwise, they 
were coded with a 3. 5, or 6, depending on the dominant type. 

This map was Incorporated Into an algorithm similar to the land-water 
algorithm. For a particular latitude and longitude the code Is retrieved 
for the box containing that location. This algorithm was used during the 
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Type (Code) 

TABLE 1 

Codes Used in Vegetation Map 
Description 

0 

Outside of Brazil Rain Forest 

1 

Humid Upper Amazon 

2 

Humid Firm Land 

3 

Humid Floodlands 

4 

Humid Plains 

5 

Mountainous Region 

6 

Savannah 

7 

Combination Typts 1,4 

a 

2,4 

9 

1 , 2 , 4 

10 

1,2 

11 

2,6 

12 

2,4,6 

13 

1 , 2 , 3 , 6 

14 

1,6 

15 

2 , 3,6 

16 

2,3 

M 

1,2,3 

18 

1,3 

19 

2,5 

20 

2 , 5,6 




SASS aiialysis to remove all data from vegetation types 3,^,5, and 6 . A 
similar approach Is recommended for a NOSS-type of system. 


4.0 CHARACTERIZATION OF g° SIGNATURE USING SEASAT SCATTEROHETER DATA 

4. 1 The SEASAT-A Scattercmeter System (SASS) 

The SASS operated using four dual-polarized (V and H) 0.5“ x 25“ fan- 
beam antennas. The antennas were pointed 1*5“ from the satellite subtrack 
to give a "star-like" illumination pattern on the surface [ 3 ] (see Figure 3 ), 
This configuration provided for pairs of nearly orthogonal measurements of 
the earth's surface. The SASS instrument transmitted a 100-watt l4.6-GHz 
signal. The backscattered signal was spread out In frequency due to the 
Doppler effect. Using 15 parallel channels (each consisting of a Doppler 
filter, a square-^law detector, and a gated integrator), the reflected power 
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FIGURE 3 


SEASAT Scatterometer Swath Coverage 
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from one antenna was sampled 61 times during a l,89“sec measurement period. 

The mean values of the 61 Integrated voltage levels of signal plus noise 
and noise alone from each of the 15 channels were transmitted to the ground. 
The power reflected from the earth's surface was determined from these vol- 
tage pairs, the calibrated gain of the receiver, and the known Integration 
time, as indicated In Figure ^4. 

Figure 5 shows an enlarged view of one of the 15 Doppler cells synthe- 
sized along the beam. As Illustrated, the Doppler-cell area Is determined 
by the 0.5* 3“dB beamwidth In the narrow-beam plane and the Doppler filter 
along the beam. Due to the 1.89-sec measurement time the Doppler cell Is 
smeared In the alongtrack direction to give the resultant cell shown In the 
lower Inset. The integrated Doppler cells are between 16-20 km wide and 
50-70 km long. 

The frequencies selected for the Doppler filters result In three near- 
nadlr resolution cells (0“ - 13® incidence angle) and twelve off-nadlr cells 
(20® - 65 ® Incidence angle). Only the off-nadir measurements were used in 
this analysis since the near-nadir returns were outside the main antenna beam. 

The value of can be determined by solving the radar equation In the 

following form: ^ 

G 2 f ( 7 ^) „ 

^ (i»7T)3 rC 

A 

where : 


= received power 
P.^. = transmitted power 
X « transmitted wavelength 


I 



H 


Lg = miscellaneous losses 
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FIGURE 5: Doppler Cell Area Determined by Antenna Pattern and Filter Response 
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Gq - peak antenna gain 

Q 

7 — ■ relative gain at center of cell 

h(f) “ freqt'cjcy response of Doppler filter 
R ■ range to call center 


Equation (l) may be solved for a as 


( 2 ) 


assuming: 


h(f) “ i 


I f„ < f < f 
0 otherwise 


where: 

A = L R<f)^, vdicro L Is the distance on earth's surface from fjj^ to f^ 

f^ ■ center frequency 

B ■ noise bandwidth (Doppler filter) 
n 

f-,f * upper and lower Ooppler-f 1 Iter cutoff frequencies 

U 

<p^ • narrow dimension antenna beamwldth (radians) 

B 

B 

f - f + -T?- 
u o 2 


A particular 0 ° Is calculated by evaluating equation ( 2 ) using measured 
Values for and P^, values for 6 q and -g— obtained from a lookup table, 
and computed values for R and A [AJ . 


4.2 Target Stability of the Anazon Rain Forest 

Initially the SASS data were processed by scientists at NASA Langley 
Research Center to produce time histories of over Brazil [4]. Figures 
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6 and 7 *vhow sample time histories of Beam I vertical-polarization measure- 
ments made during an orbit Identified as Rev 952. The radar-cross-section 
stability over the Amazon rain forest was similar to that observed during 
the SKYLAB mission. When the rain rorest results are compared with the 
observations over the ocean, three striking differences are noted; (1) o° 

Is consistently higher over the rain forest (except at the smallest Incidence 

angle), (2) 0° Is much more stable over the rain forest, and (3) the varia- 
tion In 0*^ with Incidence angle Is much less. As can be seen In both Figures 
6 and 7. SASS cells which Included part of the Amazon or one of Its major 
tributaries were easily Identified by the sudden 1-3 dB drop In 0°. 

Since the goal of the research was to determine the Amazon rain forest's 
suitability as a standard target, tests were made for both temporal and 
regional stability. Comparing data taken at the same time of day over the 
same flight line throughout the mission provided a useful test for seasonal 
effects. Figure 8 shows the results of one of these comparisons. Approxi- 
mately 5 measurements were averaged at each Incidence angle to estimate the 
mean 0°. The deviation of the mean 0° Is negligible near the peak antenna 
gain (9j * except for the September l8 measurements at ^0". This dis- 

crepancy may be due to atmospheric attenuation. At the extremes of the Inci- 
dence-angle range the deviation Is less than 0.5 dB. However, due to the 
short lifetime of SEASAT (99 days, July-October) , further research Is needed 
to determine If any seasonal effects exist. The values of the mean o° for 
a particular time of day were compared over numerous flight lines to test 
for regional effects. No consistent discrepancies were observed over the 
rain forest. As a result of these tests, the entire SEASAT rain-forest data 
set was analyzed as a unit. 
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FIGURE 8: Target Stability During SEASAT Mission 
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^.3 Incidence Angle Dependence and Diurnal Variation 

The SASS made measurements over the Amazon rain forest twice dally. 
During the early part of the mission, the Brazil crossings were between 
0500 and 0630 (local time at center of area) and between I 63 O and I 63 O, 
During the last month of the SEASAT mission, the crossings took place 
between 0900 and 1200 and between 2100 and 2li00. For the subsequent analy- 
ses the data were binned by polarization, beam. Incidence angle and time 
of day. Mean values and standard deviations of were computed for each 
of these bins. 

Plots and regressions of mean 0 ° versus Incidence angle were made 
for each beam to describe the small Incidence-angle dependence observed 
In the time histories. The sample sizes for estimation of the means varied 
approximately from 10 to 100 measurements. Figure 9 shows the results for 
Beam vertical polarization. The late-morning and nighttime data agree 
well. The results from Beams 1 and 2 Indicate that the 1 ate-af ternoon 
measurements were also consistent with the late-morning and nighttime data. 
On the other hand, the early-morning measurements were from 0.5 to 1.0 dB 
higher than data at the other times of day. 

As can be clearly seen from Figure 9, the basic angular trend Is 
adequately modeled by the equation 

a^d8) " a0 + b (3) 

where 6 Is the Incidence angle. Typical values for the Incidence angle 
coefficient a, as determined from regression analysis, ranged from 0,077 
to 0.159 dB/deg depending on beam and time of day. 
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where; 


K 



10 




Table 2 gives a complete summary of the values of a, b, K, 0^ determined by 

"it 

regression using (3) * From the summary In Table 2 some obvious differences 
are observed between beams. The data appear to be most consistent in the 
incidence angle range of 30® - 50®. The four SASS antennas had been cross- 
calibrated using the eleven Brazil passes available to NASA Langley scien- 
tists soon after the SEASAT mission [^] . Good agreement was achieved in the 
mid-angular range. The remaining differences highlight the need for the more 
sophisticated cross-calibration algorithm (using all of the SEASAT passes) 
developed here. 

The value of estimated from the regression fit agrees well with the 
reported SKYLAB results. The SASS, hoviever, shows much less deviation about 
the mean. Bars have been drawn on Figure 9 to show the scatter in the measure- 
ments. At incidence angles near the nominal pointing angle of 44®, the 
standard deviations ranged .'I'om 0.1 to 0.2 dB. At the extremes of the inci- 
dence angle range the standard deviations ranged from 0.3 to 0.6 dB. This 

A 

Both unweighted regressions and regressions weighted by number of 
cases were tried. Weighted regression results are reported here because 
they take into account the higher density of measurements In the angular 
region of greatest measurement accuracy. In most cases the differences 
in the "a" coefficients between the weighted and unweighted regression 
were within the errors of the estimated coefficients. 
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:ers for the 
SASS Deta over Brazil 


determined from a» b. 


Number 
of cases 


0.669 

33.67 

72 

0.76i» 

27.41 

246 

0.519 

38.68 

543 

0 . k 2 B 

45.30 

551 

0./<35 

56.31 

166 

0.391 

49.00 

518 

0.663 

32.17 

442 

0.77^ 

28.09 

437 

0,874 

29.30 

103 

0.767 

27.91 

603 

0.366 

55.61 

836 

0.582 

42.29 

1 40 

0.510 

38.65 

589 

0.595 

32.95 

538 
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tncrease In the standard deviation Is probably due to a decrease In measure- 
ment precision caused by weaker signals at angles well away from the peak 
antenna gain, and also may be caused by smal'i beam-polntlng-angle errors 
that have more Importance where the antenna gain varies more rapidly with 
angle. 

The regression equations determined from the data were used to estimate 
values of 0° at 1*0'' Incidence angle which are plotted In Figure 10 as a 
function of time of day for Beam 2. The early morning return Is 0.6 - 0.9 
dB higher than at the other times of day. Table 3 summarizes these estimates 
for all of the beams. When this effect was first noted In the original 
analysis of the first 11 passes, all of the data had been screened using 
GOES Imagery to eliminate cases containing significant cloud cover. In the 
present analysis the data were not screened. However, since this was the 
dry season few of the measurements should have been corrupted. 

Initially It was believed that, because of the random orientation of 
the seatterers In the uniformly vegetated rain forest, the values of 
would be Insensitive to polarization. While most of the measurements were 
vertically polarized, enough horizontal polarization data were taken to 
allow a comparison. On Figure 11 are plots of both V-pol and H-pol data for 
Beam 1. As expected the results appear to be Insensitive to polarization. 

No noticeable azimuthal-angle variation was observed In the rain-forest 
data. This Isotropy can be demonstrated by comparing the 0900 - 1200 V-pol 
data from Beam 2 and Beam 4. While the beams differ 180° in azimuth, the 
measurements show agreement within O.l - 0.2 dB. This may be seen by over- 
laying Figures 9 and 1 1 . 
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FIGURE 10: Estimated a'' From Regression Fit 
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FIGURE 11: Comparison of VV-pol and HH-pol Data 
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TABLE 3 

Summary of Diurnal Effect at Incidence Angle 
Vertical Polarization 


Time of Day 

ESTIMATED 

SCATTERING 

COEFFICIENT (dB) 

Average 

for 

A 1 1 Beams 

Beam I 

Beam 2 

Beam 3 

Beam 4 

05G0 - 0630 

-6.91 

-6.70 

-6.3l» 

-6.1<7 

- 6.60 

0900 - 1200 

-7.A9 

-7.62 

-7.39 

-7.l»l 

-7.48 

1630 - 1830 

-7.33 

-7.^0 




-7.36 

2100 - 2 I 4 OO 

-7.53 

-7.31 

-7.49 

-7.5^ 

-7.47 


SEAS AT Summary 

The data obtained from SEASAT-1 have shown the potential of the Amazon 
rain forest for use as a standard target. Since the diurnal effects appear 
to be limited to early morning, data from other times of day could all be 
used together to build a data base for calibration. The small scatter (O.l - 
0.6 dB) In the mean values of Indicates the remarkable stability of this 
region and hence its usefulness as a standard target. A basic character- 
ization of the incidence-angle dependence can be made using a simple straight- 
line model of vs 0. 

The differences between the regression coefficients shown In Table 2 
for the different beams cannot be due to accual 0 ° differences because Beams 
1 and 2 covered the same area, as did Beams 3 and Hence these differences 
must be In antenna patterns or In pointing-angle differences. This points 
o‘jt the need for use of such a "standard target" In future multibeam scat- 


terometers. 
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5.0 DEVELOPMENT OF OFF-NADIR ANTENNA BIAS CORRECTIONS FOR NOSS 

5 . 1 Description of NOSS Scatterometer System 

The NOSS Scatterometer System (SCATT) design Is for a second generation 
fan-beam scatterometer using’ experience gained from SEASAT (Section ^.1). 
Figure 12 shows the proposed SCATT ground-track coverage. Besides the four 
beams used In SASS, two antennas are pointed at 65“ from the satellite 
orbit plane to aid In wind-direction alias removal. All of the antennas are 
dual -polarized (V and H). The SCATT Instrument transmits a lOO-watt )4.0- 
GHz signal. The backscattered signal Is Doppler filtered Into 60 resolution 
cells with each cell approximately 10 kro x 10 km. The off-nadIr bias 
correction algorithm Is to be applied to the outermost 50 Doppler cells of 
each beam. 

5.2 Development of Maximum Likelihood (ML) Algorithms 

to Estimate Off-Nadir Bias Parameters 

During SEASAT analysis, Brazilian rain-forest data were first used to 
correct for antenna biases In the off-nadIr cells. Estimates for relative 
bias and pointing angle were determined by comparing data from a limited 
number of passes. For NOSS a formalized algorithm was desired to allow 
automation of the estimation of these parameters by using all available 
SCATT rain-forest data. The next subsection (5.2.1) describes In detail 
the derivation of a maximum likelihood (ML) approach to estimate relative 
bias and pointing angle. Once pointing angle Is determined, a ML technique 
may be used periodically to re-estlmate relative bias. This algorithm Is 
described in Section 5*2.2. 


*Thls ML approach Is an adaptation of the approach used in the SEASAT 
V/lnd Vector Algorithm developed by F.J. V/entz (Remote Sensing Systems, Inc.). 
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FIGURE 12: Proposed NOSS Footprint 



28- 


Original PAGE fS 
OF POOR QUALITY 


5 > 2 .] Relative Bias and Pointing Angle Estimation 

For an operational system the radar parameters, Pj and Gq In equations 
(1) and (2) are carefully specified In the design. The relative antenna 

Q 

gain, Is determined from pre-launch antenna pattern measurements and the 
specified pointing angle, 0 . However, as was observed v/lth the SASS, the 
actual values for P^, and 0p after launch can differ from the specifi- 
cations [4], 

Effects of these differences can be seen by writing equation (1) for 
both the designed and actual specifications of P.j., Gg, 0p and the corres- 
ponding a° Inferred from equation (2). 


P« “ 
R 




r"' 


where the subscript D refers to the designed values. 

P U^G Q^[G/Gn(0 „)] a° A 
„ _ A^ A° ° aP a 


R 




where the subscript A refers to actual values. 


Equating the above expressions and solving for Og In terms of 


yields 


P[G/G„ 


:g/o„ (e^p)r 


(5) 


A ^A^ 

where a = relative bias = : 

P T 0 

d' d 


By substituting a standard reference target, Og, for the value of a^, 
the relationship In equation (5) can be used to estimate the relative bias, 


a, and the actual pointing angle 0 


Since the measurements, Op., are noisy, 
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maximum-! ike 1 1 hood estimation is used to determine <x and 9 . 

aP 

Formulation of the maximum-] ike) i hood estimation may be given as: a 

vector of noisy observations, z, exists which depends probabilistically on 
a parameter of Interest, B. The M.L. estimate of B is then given by the 
condition 



■ maximum 



iz) 


( 6 ) 


where: 


f is the Joint conditional probability density function. 
5/B 

A 

bj^^(z) is the maximum-likelihood estimate of B. 


This approach can be applied to estimate 
assumptions : 


a and 9 with the following 


(1) Radar measurements are independent and may be described 
as a Gaussian process 

(2) The noisy observations, z, are the parameter 

(3) The parameter of interest, B, is given by the relation: 


a [G/G (0 )]^ 

a P 

[g/Gq (e^]2 


B(a 0 p) 
A 


(4) 


G/GoOp) 


is known from pre-launch measurements. 


Applying the above conditions, the conditional probability in expression (6) 
may be re-written as: 
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f. [{a° ; £ - 1,N}/D(a,e J] 

i P, A U . □ 


a°/B 


whe re : 


N 

n e 
A" I 


-[0°P - B(a,6 p)]^/2(AaJ)2 
A 


N y — 

n /2TT (Aa°) ^ 


( 7 ) 


Since the purpose is to select values of a, 6 ^ that maximize expression 
(7), two simplifications can be made. First, the quantity, n, may be Ignorec 
since It Is not a function of the parameters a,0 . Second, the logarithm 

aP 

of the likelihood function may be maximized, since the same values a»6 

r 

will maximize both the likelihood function and the log-1 Ikel Ihood function. 
Therefore, the condition to satisfy Is: 


-1/2 ? - B)' 

A=1 ^ 


max 


( 8 ) 


where B » B(a,0 ). 

aP 

Since the function G/GQ(0p) Is stored as a look-up table and very 
difficult to express analytically, the log-1 Ikel ihood expression In (8) Is 
maximized In two steps. First a migrating search Is performed to approxi- 
mately center a 3 X 3 matrix about the solution. Then the maximum of a 
bivariate Interpolating polynomial fit to points on this matrix Is found 

A A 

to obtain a,0 . 

P 

The search Is done by defining a log- 1 Ikel Ihood matrix, j, as follows 




( 9 ) 
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where: 

+ lAa; I ■ -1 ,0, 1 
®Jp “ ®op + J 


Since the actual values for and Gq should be reasonably close to 

the design specifications, Initially Is chosen to be I.O, Similarly, 

the design-specified pointing angle, 6 , Is chosen for the Initial yalue 

D° 

of 0Qp. Reasonable values for Aa,A0p appear to be 0,2 and 1,0 respectively. 
The values for gj^j are then compared to find the maximum. If the center 
point, pQ Q, Is not the maximum, a new matrix Is formed as follows: 


“I -“nax + ' 

®ID " ®D J ^ 

^ ^ ^max ^ 

where: 


a ^,0 are the a,,0,_ corresponding to the maximum of the previous 
max' p I j P r = r 

i»iaA 

matrix. 

Searching continues until a matrix Is formed which has the maximum at the 
center. 

A A 


Estimates for relative bias, a, and actual pointing angle, 0 , are found 

P 

by fitting a bf-varlate, quadratic Interpolating polynomial to the values 

A A 

of g, Values of a and 0 which maximize the polynomial are then deter- 
• tJ P 

mined analytically. The Interpolating polynomial Is of the form: 


g(ci,0p) 


+ bg + cy^ dy + e0y + f 


( 10 ) 


where; 
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a 




c 


^ 0,-1 

2 



d 




® " ^ 0,0 " ^ 1 ,0 " ^ 0,1 ,) 

^ “ 90,0 


where g, , are the values of the likelihood matrix. 

I 

The maximum Is found by taking the partial derivatives of the Inter- 
polating polynomial and setting them equal to zero. When the resultant 

A A 


equations are solved for a and 0 we get 

r 


a, 


ML 


a. 


, Act(ed - 2bc) 

T 5 ' 

Ijac - e 


Ola) 


A 


e 

pml 



A0 (be - 2ad) 
kac - e^ 


(lib) 


where a,b,c,d,e are defined as above. 


5.2.2 Long-Term Monitoring of Relative Bias 

Periodic estimates of the relative bias during the mission will provide 
a means to correct for long-term drifts In transmitter power. If the 

A 

pointing angle Is not expected to change once 0^ Is determined, a simpli- 
fied version of the algorithm In Section 5.2.1 may be used to estimate 


relative bias. 
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In this case the parameter of interest contains only one unkno^^Ot 
the relotlvo bias a. The same procedure followed In Section 5*2, 1 is used 
to solve this one-parameter case with a likelihood vector^ gp replacing 
the motrix g, The vector Is defined as: 


-4 I 

^ !l-l 




Oi2. 




i - -1,0,1 


( 12 ) 


where: 


.0 


Ot, 




■ the NOSS SCATT measurement 
** the standard torget 
« Qq + lAct 

•• previously determined relative bias 


Vector gj Is then fitted with a quadratic Interpolating polynomial 
of the form: 


g(a) 



+ s8 + t 


(13) 


where : 

6 

r 

s 


An 




t - go 

g| the values of the vector calculated in equation (12). 

By differentiating equation (13) and setting the result equal to zero, 

A 

the maximum likelihood value of a Is found as; 
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5.3 Implementation of Off~Nad[r Bias Algorithm for NOSS SCATT 

The flow chart in Figure 13 shows the relationship between the off- 
nadir bias algorithm and SCATT data processing. Earth-located a® data 
are screened to locate data over the Amazon rain forest. Then the Amazon 
map described In Section 3.0 Is used to remove measurements that Include 
either the Amazon river or one of Its principal tributaries. Remaining 
rain-forest data are then supplied to the off-nadIr bias algorithm. 

Actual estimates of the relative bias and pointing angle are performed 
off-line from the main SCATT processing. This Is necessary for two reasons: 
first, Brazil passes for several days are needed to create a sufficient 
data base for estimation; second, considerable human Interaction and evalu- 
ation is needed to Interpret the results. 

FlgiJre 11, shows a flow chart of the processing required in the off- 
nadlr bias algorithm. Inputs to the algorithms are the SCATT data from 
the rain forest and the corresponding brightness temperatures Tg measured 
by the Large Antenna Multichannel Microwave Radiometer (LAMMR). Co-located 
LAMMR data are used to flag and remove potential rain-corrupted SCATT 
measurements. A data base Is then created for use with the maximum- 1 I kel I- 
hood procedure described In Section 5.2. Estimates for relative bias and 
antenna pointing angle are computed for every combination of beam, cell and 
polarization. Approximate error bounds are determined using the results of 
several estimations. The following subsections describe in more detail the 
processing required. 

5 . 3.1 Co-Location of SCATT and LAMMR 

To use LAMMR 37 GHz data to flag rain-corrupted SCATT measurements, 

To measurements must be co-located with individual SCATT 0 ° measurements. 
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FIGURE 


3 


Proposed Processing of NOSS Scatterometer (SCATT) Data to 
Obtain o° 
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SC ATT DATA LAMMR DATA 



Relative Biases, Pointing Angles 


FIGURE U: Processing Required for Off-Nadir Bias Correction 
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The Input data streams for both SCATT and LAMMR are assumed to be gridded 
Into 50 km blocks and assigned a set of coordinates (I,j). Co-locatlon 
Is accomplished by: first, sorting both LAMMR Tg's and SCATT o°'s Indepen- 

dently by grid coordinates; second, for each coordinate pair (l,J) contain- 
ing SCATT data the corresponding LAMMR grid point Is located. Within the 
50 km grid each SCATT cell Is matched with the closest LAMMR cell. 

5.3«2 Flagging Rain Cells 

Co-located LAMMR 37 GHz Tg's are used to flag potential rain-corrupted 
SCATT cells. The brightness temperature of the Brazilian rain forest Is 
very high due to high emisslvlty and surface temperature. However, In the 
presence of rain only the top portion of the rain cloud and not the 
forest is seen. Hence It is believed that the temperature measured should 
be lower when rain occurs than In clear sky conditions. Rain flags will be 
set whenever the measured temperature falls below a fixed cut-off temper- 
ature, T^^^. 

Figure 15 shows the form of the 37 GHz brightness temperature histogram 
expected from the Amazon rein forest. The lobe at the lower brightness 
temperature corresponds to measurements made In rain conditions, while 
the lobe at higher brightness temperature corresponds to clear-sky conditions. 
A preliminary value for T^^^^ i,.-''' possible be made from SEASAT radiometer 
measurements. However, ft will be Important to make a number of histo- 
grams with LAMMR data to confirm the assumed distribution shown In Flgiire 
15 and to properly set T^^^. 
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5*3*3 Data Base Creation 

Due to the largo volume of data produced by tho SCATT, estimation of 
relative bios and pointing angle Is performed by using moan values com- 
puted for each pass over Brazil. Indlvlduol clear-sky SCATT measurements 
are grouped by beam, coll number and polarization. At the end of each Brazil 
pass the mean value, standard deviation, and average Incidence angle is 
computed for each group containing more than twenty measurements. These 
estimates are also binned by beam, cell number and polarization. Estimation 
of relotivo bias and pointing angle takes place when each mean value bin 
has data from at least 10 passes. Since the bins will fill at different 
rates depending on the location of the satellite subtrack, provision Is made 
to store up to 20 values In each. 

5 . 3 .^ Standard Target Creation 

Tho maximum likelihood procedure described in Section 5.2 requires a 
standard target, cJg(O), defined over the entire off-nadIr Incidence angle 
range. Determination of Og will require consistent human evaluation and 
deci s Ion. 

Initially plots and regressions of the mean value of 0 ° versus Incidence 
angle will be made from the data base der.crlbed above. These plots and 
regressions will be compared with SEASAT results described In Section l<.0 
and Appendix B. If there Is good comporlson In trend, the SEASAT model, 

Ojg ■ aO + b, will be adopted with an appropriate level shift. If the trend 
of the SCATT data Is significantly different, a new model determined from 
the analysis of several passes may be used. The final procedure used In 
creating the standard target will be determined once the Initial comparisons 
are made. 
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5.3,5 Relative Bias and Pointing Angle Estimation 

Early In a future mission with the NOSS design will be the Implementation 
of the HL technique to determine both relative bias and pointing angle 
described In Section 5>2,2, Once a suitable valu'i for pointing angle Is 
determined for each beam, the simplified algorithm described In Section 5.2.3 
may be used to monitor long-range transmitter drift. Both of these algorithms 
can be Implemented using the same Inputs so switching between algorithms 
should be trivial . 

In this Implementation the noisy data, in equations (8) and (12) 

r 

are the mean values stored In the data base (Section 5>3.3). Separate 
estimates will be made for each combination of beam, cell and polarization. 
Repeated estimates should be made* using subsequent data. Final estimates 
and bounds may be calculated by computing the means and standard deviations. 


6.0 CONCLUSIONS AND FUTURE RESEARCH 

Detailed analysis of SEASAT scatterometer data has confirmed the suita- 
bility of the Amazon rain forest for use as a standard calibration target. 

A simple straight line model for Incidence angle variation appears to be 
quite adequate. Useful approximations for the model parameters may be deter- 
mined from the SASS data. The accuracy of these approximations Is limited by 
remaining biases between SASS antennas and the lack of Independent rain-forest 
measurements for absolute calibration. 

Max I mum- 1 Ike) 1 hood estimation algorithms have been developed to correct 
off-nadir antenna biases for spaceborne scatterometers by utilizing raln~ 
forest data. Application of those algorithms for the proposed NOSS scatter- 
ometer was straight-forward. Presumably these techniques could be adapted 
easily for calibrating other future systems. 




Three Important areas of research need to be undertaken to extend the 
usefulness of Amazon rain forest data to correct off-nadIr antenna biases. 
First, measurements mur.t be made throughout the year to determine the exis- 
tence of any seasonal effects. The SEASAT measurements were made during 
only the dry season. Second, research needs to be done to confirm the pro- 
posed algorithm to screen out data corrupted by thick clouds and rain. The 
multifrequency radiometer measurements made with the NIMBUS 7 satellite 
could potentially contribute valuable information In both of these areas. 
Third, a calibration program using aircraft overflights Is needed to better 
determine the absolute value of . 
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CBRZLND lAND/WATER ROUTINE FOR BRZILIAN RAIN FOREST 

0 

SUBROUTINE LDBRZL 
PURPOSE 

THIS ROUTINE LOADS THE LAND WATERMAP 

INTEGER ARRAY ( 6 , 1 0 , 64 ) , CODE 
REAL GLAT(6)/2. ,0. ,-2. ,-4* ,-6. ,-8./ 

REAL GL0N(l0)/74. ,72. ,70. ,68. ,66. ,64* , 

& 62. ,60, ,58. ,56./ 

REAL LAT.LONG 

READ LAND/WATER MAP 

DO 10 1=1 ,6 
DO 10 J=1 ,10 

READ (01 ,900) (ARRAY(I,J,K) ,K=1 ,64) 

900 FORMAT (4X, 818) 

1 0 CONTINUE 


RETURN 


ENTRY LNDH20(TLAT,TL0NG, 

CODE) 

PURPOSE 


THIS ROUTINE FETCHES 
CORRESPONDING TO AN 
LONGITUDE 

A LAND/WATER CODE 
INPUT LATITUDE AND 

INPUT ARGUMENTS 


TLAT * 

TLONG * 

LATITUDE (DEG) 

EAST LONGITUDE (DEG) 

OUTPUT ARGUMENT 


CODE * 

LAT = TLAT 
LONG = TLONG 

LAND/WATER CODE 

LONG = 560.0 - LONG 



ILAT = IFIX( (2.-LAT)/2. ) + 1 

JLON = IFIX (74.-L0NG)/2. ) + 1 

ITEP = IPIXC (GLAT(ILAT) - LAT)/0.25 ) 

JTEP = IPIX( (GLON(JLON) - LONG) /O. 25) + 1 

KELE = ITEP*8 + JTEP 

CODE = 1 + ARRAY{ ILAT, JLON, KELE) 

RETURN 

END 

CBRZVEG VEGETATION CODE ROUTINE FOR BRAZILAIN AREA 

C 


0000001 0 
00000020 
00000050 
00000040 
00000050 
00000060 
00000070 
00000080 
00000090 
00000100 
00000110 
00000120 
00000150 
000001 40 
00000150 
000001 60 
00000170 
000001 80 
000001 90 
00000200 
00000210 
00000220 
00000250 
00000240 
00000250 
00000260 
00000270 
00000280 
00000290 
00000500 
00000510 
00000520 
00000550 
00000540 
00000550 
00000560 
00000570 
00000580 
00000590 
00000400 
00000410 
00000420 
00000450 - 
00000440 
00000450 
00000460 
00000470 
00000480 
00000490 
00000500 
0000051 0 
00000520 
00000550 
00000540 
00000550 
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SUBROUTINE LDVEO 
PURPOSE 
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THIS ROUTINE LOADS THE VEGETATION MAP 
INTEGER VEGTBL( 5 ,3,64), VGCODE 
REAL LAT,LNG 

REAL VLAT{3)/2.0,-2.0,-6.0/ 

REAL VLNG(5)/74. 0,70. 0,66. 0,62. 0,58.0/ 

READ VEGETATION MAP 


DO 10 I = 1 ,5 
DO 10 J s: 1 ,3 

READ(02,900)(VEGTBL{I,J,K) ,K=1 ,64; 
1 0 CONTINUE 
900 FORM AT (4X, 81 5) 


VEGCD 


RETURN 

COMPUTE CORRECT VEGETATION CODE 


ENTRY VEGCD( TLAT,TLNG, VGCODE ) 


PURPOSE 


THIS ROUTINE INPUTS A LATITUDE AND A 
LONGITUDE. WITH THESE A VEGETAION 
CODE IS DETERMINED 

INPUT ARGUMENTS 

TLAT * LATITUDE (DEG) 

TLONG * EAST LONGITUDE(DEG) 

OUTPUT ARGUMENT 

VGCODE * VEGETATION CODE 

LAT = TLAT 
LNG = TLNG 

CONVERT TO WEST LONGITUDE 

LNG = 360.0 - LNG 

DETERMINE VEGETATION CODE 

ILAT = IPIX( (2.0 - LAT) / 4.0) + 1 
ILNG = IPIX( (74.0 - LNG) / 4.0) + 1 

INDEXX = IPIX( (VLAT(ILAT) - LAT) / 0.5) 

INDEXY = IPIX( (vLNG(iLNG) ~ LNG) / 0.5) + 1 


00000560 
00000570 
00000580 
00000590 
00000600 
00000610 
00000620 
00000630 
00000640 
00000650 
00000660 
00000670 
00000680 
00000690 
00000700 
00000710 
00000720 
00000730 
00000740 
00000750 
00000760 
00000770 
00000780 
00000790 
00000800 
0000081 0 
00000820 
00000830 
00000840 
00000850 
00000860 
00000870 
00000880 
00000890 
00000900 
0000091 0 
00000920 
00000930 
00000940 
00000950 
00000960 
00000970 
00000980 
00000990 
00001 000 
00001010 
00001020 
00001030 
00001040 
00001050 
00001060 
00001070 
00001 080 
00001090 
00001 100 
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INDEXZ * INDEXX 8 + INDEXY 000011 1C 

0 00001 1 2C 

VGCODE s VEGTBL( ILNG , XLAT , INDEXZ ) 00001 1 5C 

C 00001 14C 

RETURN 00001 1 5C 

END 00001160 
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